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Abstract 

The random fiber bundle (RFB) model, with the strength of the fibers distributed uniformly 
within a finite interval, is studied under the assumption of global load sharing among all un- 
broken fibers of the bundle. At any fixed value of the applied stress a (load per fiber initially 
present in the bundle), the fraction Ut{(j) of fibers that remain unbroken at successive time 
steps t is shown to follow simple recurrence relations. The model is found to have stable 
fixed point f/*(cr) for applied stress in the range < a < cTc, beyond which total failure of 
the bundle takes place discontinuously (abruptly from U*{ac) to 0) . The dynamic critical 
behavior near this has been studied for this model analysing the recurrence relations. 
We also investigated the finite size scaling behavior near ctc. At the critical point a = ctc, 
one finds strict power law decay (with time t) of the fraction of unbroken fibers Ut{<Jc) (as 
t — > cxd). The avalanche size distribution for this mean-field dynamics of failure at cr < Uc has 
been studied. The elastic response of the RFB model has also been studied analytically for a 
specific probability distribution of fiber strengths, where the bundle shows plastic behavior 
before complete failure, following an initial linear response. 
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I. Introduction 

A typical relaxational dynamics has been observed in a strained random fiber bundle (RFB) 
model [1-8] where N fibers are connected in parallel to each other and clamped at their 
two ends and the strength of the individual fibers has some particular distribution (white, 
Gaussian or otherwise). In the global load-sharing approximation ||l], ||, at any instant, the 
surviving fibers all share equally the external load (irrespective of their proximity etc. of 
the fiber to failed fibers etc.). Initially, after the load F is applied on the bundle, a fraction 
of the fibers having strength less than the applied stress a = F/N fail immediately. After 
this, the total load on the bundle redistributes globally as the stress is transferred from 
broken fibers to the remaining unbroken ones. This redistribution causes secondary failures 
which in general causes further failures and so on. After some relaxation time r, which 
depends on cr, the system ultimately becomes stable if the applied stress a is less than or 
equal to a critical value cTc, and beyond which (a > cTc) all the fibers break and the bundle 
fails completely. Although the local load sharing might be more realistic, we study here the 
global load sharing model because of its simplicity. The study of the scaling properties of the 
dynamics of the fiber bundle model systems is expected to be extremely useful in analysing 
the statistics of fracture and breakdown in real materials, including earthquakes 0, p!0|| . 

In this paper, we report on the critical dynamics of the RFB model in the global load- 
sharing case, assuming uniform distribution of threshold strength of the fibers (up to a 



cutoff), in particular at the critical point cTc. In a previous paper |[Tl|, we have solved 
the dynamics of the model, showing a novel critical behavior as the stress a approaches a 
certain value ac, we had derived there the expressions for the breakdown susceptibility x 
and the relaxation time r under a stress cr < Uc and showed that both the quantities diverge 
following power laws as a approaches cTc from below. Here we define an order parameter for 
the transition from a state of partial failure of the bundle to a state of total failure and also 
show that at the critical stress dc, the dynamics follows a precise and strict power law. From 
the finite size dependence of cxc and the order parameter we have identified the correlation 
length exponent of the system. We have studied the avalanche size statistics in the model 
as well. Considering a modified (uniform but shifted from the origin) distribution of fiber 
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strengths we have studied analytically the elastic-plastic deformation characteristics of 
the RFB model. 

II. The model 

The RFB model consists of elastic fibers clamped at two ends (Fig. 1), where the failure 
stress of the individual fibers are distributed randomly and uniformly in the interval between 
and 1 (white or uniform distribution; Fig. 2). Global load sharing is assumed; i.e., the 
applied load on the bundle is equally shared among all the existing intact fibers. This 
assumption neglects 'local' fiuctuations in stress (and its redistribution) and renders the 



Fig. 1: The RFB model consists of N fibers. The bundle is subjected to a load F. Assuming linear elasticity, 
with identical elastic constant k for each of the fiber up to the breaking, the load F can be expressed as NkS 
where 6 denotes the strain for the fibers until any of them breaks. The breaking strengths of the fibers are 
assumed to be random, as discussed later. 

p(o) 

1 1 



Fig. 2: The simplest model considered here assumes uniformly random distribution or white distribution 
p(cr) for the strength of the fibers up to a (normaHsed) cutoff strength. This distribution gives the recurrence 
relation (1). 



model as a mean-field one. We work with the fraction Ut = Nt/N; Nt being the number of 
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fibers remaining intact after t time-steps and A^t=o = N. 

With the application of any small load F {— a N, with o" ^ 1) on the bundle, an initial 
stress a (load per fiber) sets in. At this first step therefore, aN number of fibers break, leaving 
NUi{a) = N{1 — a) unbroken fibers. After this, the applied force is redistributed uniformly 
among the remaining intact fibers and the redistributed stress becomes F/[NUi{a)] — a/{l — 
a). Some more fibers, for which the strengths are below the value of the redistributed stress, 
fail thus leaving NU2{(t) — N[l — a/{l — a)] unbroken fibers. This in turn increases the 
redistributed stress and induces further failures. Consequently, as the stress per fiber cr(t) 
at time t is given by F/NUt = a/Ut and the surviving fraction is given by 1 — a/Ut (see Fig. 
2), Ut{a) follows a simple recurrence relation 

^«<'') = ^-ci^- « 

III. Breaking dynamics of the RFB model 

The recurrence relation (1) has the form of an iterative map Ut+i — Y{Ut). Its fixed point 
U* is defined by the relation: U* = Y{U*) and from eqn. (1) one gets 

U%a) = l±{a,-ay/';a^ = \. (2) 

The quantity U* must be real valued as it has a physical meaning: it is the fraction of the 
original bundle that remains intact under a fixed applied stress a when the applied stress 
lies in the range < cr < Uc- For a > ac the map does not have a real- valued fixed point 
and as can be seen from (1), the dynamics never stops until Ut — when the bundle breaks 
completely. Since it requires that \dY/dU\u*{a) <1 for a fixed point U*{a) to be stable, for 
each value of a the value of U* with the positive sign in eqn. (2) represents a stable fixed 
point (or attractor) while the value of U* with the negative sign in eqn. (2) represents an 
unstable fixed point (or repeller). It may be noted that the quantity U* — 1/2 behaves like 
an order parameter that determines a transition from a state of partial failure (cr < dc) to a 
state of total failure (cr > (7c) : 

O = U* -1/2^ {a,- af;(3^^. (3) 
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IV. Critical behavior 

(a) For cr < (Tc 

To study the dynamics away from criticality (a —>■ ac from below) , we replace the recurrence 
relation (1) by a differential equation 

dU U"^ -U + a 



dt 



U 



(4) 



Close to the fixed point we write Ut{cr) = U*{a) +e which, following eqn. (4), gives jll 



Ut{a)-U%a)^exp{-t/T), 



(5) 



where r 



|(crc — a) + 1 . Approaching ac from below we get 
r oc ((Tc - a) 



1 

"=2 



(6) 



as the relaxation time of the model and it is found to diverge following a power-law as a —>■ 
from below. Although, we have used here the continuum-time version (4) of the recurrence 
relation to evaluate the relaxation time (r), we have checked numerically as well from the 
discrete-time recurrence relation (1) and obtained the same exponent value. 

One can also consider the breakdown susceptibility x, defined as the number (fraction) 



of fibers that break due to an infinitesimal increment of the applied stress [11 

dU*(a 



X 



da 



[CTr 



a) 



(7) 



from equation (2). Hence x too diverges as the applied stress a approaches the critical value 
ac = \. Such a divergence in x had already been observed in the numerical measurements 

(b) At a = ac 

At a = ac the fraction of fibers surviving is U*{ac) = \ and \dY / dU\u*{c7a) = 1 which suggests 
that the system will take infinite time to reach the fixed point at ac- From the recurrence 
relation (1) it can be shown that this decay of the fraction Ut{ac) of unbroken fibers that 
remain intact at time t follows a simple power-law: 
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Ut = -(1 



1 



2'" ■ t + 1' 
4 



(8) 



starting from Uq — 1. For large t {t ^ oo), this reduces to Ut — 1/2 oc t a simple but 
strict power law. 

V. Finite size effects and correlation length exponent 

For a finite bundle of N fibers, the recurrence relation (1) will be replaced by 

ivV 



Nt+i{a)^N 



N, 



(9) 



where [xj denotes the greatest integer less than or equal to x. Here the fixed point is 
obtained when A^t+i — Nf — N* and the value of N* is bounded by the relation 



N 



1/2 



<N*<-{N + 1) + 



.1/2 



(10) 



which clearly depends on the finite size of the system. Consequently the effective critical 
point (Tc{N) for the finite RFB model is bounded as: 

1 .... 1. 1 



-<a.(Ar)<-[i + -f. 



(11) 



It follows from eqn. (10) that, at ctc = |, the fixed point value for a finite bundle decays 
with the initial bundle size following a power law 



(12) 



Since the quantity f/* — 1/2 in eqn. (3) behaves like an order parameter for a phase transition, 
the corresponding quantity in a finite bundle of N fibers would be 



2 =U\ 



1 



N 2- 

Expressing the correlation length as ^ oc {ac — a)^'^ in the infinite system and combining it 
with eqn. (3) for a finite size system (where ^ ^ N), the finite size scaling behavior can be 
written as 

uu^c) -l-^ N-^/r (14) 

Since /3 = 1/2, as obtained earlier from eqn. (3), we get v = Ihy comparing eqn. (14) with 
eqn. (12). 



(13) 
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VI. Avalanche size distribution 

We now study the avalanche size distribution in this mean-field model. If one considers 
strictly uniform strength distribution of the fibers in this model, one can not meaningfully 
approach the failure point by breaking the weakest fiber and looking for the avalanches of 
successive failures of the fibers, following the avalanche definition of Hemmer et al |^, ||. If 
we apply this definition in the above (restricted) model, we will end up with only two distinct 
sizes of avalanches: {N/2) avalanches of unit size and one avalanche of size {N/2). This will 
occur due to the perfectly uniform strength distribution of the fibers (with the successive 
strength of fibers differing by 1/A^). To work therefore with a more general definition of 
avalanche , we increase the external load on the bundle steadily such that the external 
load F increases by an equal amount {dF = Nda) at each step (cf. 0). This ensures 
the bimodal, yet decreasing, distribution function mentioned above to become a smooth 
(decaying) function. Operationally also, this procedure is quite common and can be applied 
to different cases and to bundles with different types of strength distribution p{a) of fibers. 
Here, the fraction of fibers m which eventually fail due to this increase in load or stress may 
be considered as the avalanche size: 

m='^-M = l-W{a). (15) 
da 

With U*{a) from (2) we get 

ac — o" ~ m~^. (16) 

If we now define the avalanche size probability distribution by P(m), then P{m)Am measures 
Ao", the number of times one has to change a (by da) to get a change Am along the m versus 
a curve in (16). In other words, 

p(m) = 4^ ~ m"''; r/ = 3. (17) 
dm 

This mean-field result for P{m) (power law decay with exponent i] = 3) is obtained 
here for global load sharing and uniform fiber strength distribution when the external load 
is increased by a fixed amount. We have checked this result numerically for different da 
values ( = 1/A^) for bundles with 50,000 fibers having both strictly uniform and uniform- 
on-average strength distributions. The results are shown in Fig. 3. The earlier numerical 
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results of Moreno et. al. for Weibull type distribution of fiber strength also confirms the 
relation (16), which implies that the cumulative distribution decreases with avalanche size 
m as m~^, in agreement with (17). 




m 

Fig. 3: The log-log plot of the average avalanche size distributions P{m) against m for N = 50,000 
with da = 1/N for strictly uniform fiber strength distribution (cross) and uniform-on-average fiber strength 
distribution (averaged over 501 bundle reaHsations; plus). The dotted fine has a slope i] = —3.0, representing 
eqn. (17). The inset shows the avalanche size distributions for uniform-on-average fiber strength distribution, 
when (a) the external load increases by a fixed amount da — 1/N (plus) and (b) the avalanches are triggered 
by breaking the next weakest fiber (star). The dotted and dashed lines in the inset correspond to 77 = —3.0 
and 77 = —2.5 respectively. 

This result (17) for the avalanche size distribution P(m) is therefore valid for other 
distributions of fiber strength (cf. 0), when the avalanche size is defined through (15). 
If one looks for the statistics of avalanches initiated by breaking the next weakest fiber in 
bundles with uniform-on-average fiber strength distribution, as in Hemmer et al then 
one gets r] = 5/2. This is shown in the inset of Fig. 3, where the avalanches are defined 
in both ways: with fixed increase in a (giving r] = 3.0) and by breaking the next weakest 
fiber (giving rj = 2.5). The difference in the above exponent values therefore originates 
from different ways of defining the avalanches; in our method of defining avalanches here. 
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the external load on the bundle increases uniformly, while in the other method the external 
load increase has intrinsic fluctuations due to the randomness of the flber strengths and the 
restriction on initiating the avalanches by breaking only the next weakest flber. 

VII. Plastic deformation and stress-strain relation 

One can now consider a slightly modifled strength distribution of such a flber bundle, showing 
plastic-deformation characteristics [||, |^. For this, we consider a RFB strength distribution, 
having a lower cutoff. Until failure of any of the flbers (due to this lower cutoff) , the bundle 
shows linear elastic behavior. As soon as the flbers start failing, the stress-strain relationship 
becomes nonlinear. The dynamic critical behavior remains essentially the same and the static 
(flxed point) behavior shows elastic-plastic deformation before rupture of the bundle. 

Here the flbers are elastic in nature having identical force constant k (see Fig. 1) and 
the random flber strengths distributed uniformly in the interval [ai, 1] with cr^ > 0; the 
normalised distribution of the threshold stress of the flbers thus has the form (see Fig. 4): 



0, < a < (JL 



(18) 



P(a) 




Fig. 4: The fiber breaking strength distribution p{(j) considered for studying elastic-plastic deformation 
behavior of the RFB model. This distribution gives the recurrence relation (19). 

For an applied stress a < none of the flbers break, though they are elongated by an 
amount 5 = cr/n. The dynamics of breaking starts when applied stress a becomes greater 
than ai- Now, for a > ctl the fraction of unbroken flbers follows a recurrence relation (for 
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p{a) as in Fig. 4): 

Ut+iia) = 1 - 
which has stable fixed points: 



NUt{a) 



1 



(7L 



(19) 



2(1 - (Tl) 



1+1- 



a 



1/2 



',crc 



1 



4(1 - <Jl) 



(20) 



The RFB model now has a critical point cTc = 1/[4(1 — ctl)] beyond which total failure of the 
bundle takes place. The above equation also requires that ctl < 1/2 (to keep the fraction 
U* < 1). As one can easily see, the dynamics of Ut{a) for a < ac and also at a — ac remains 
the same as discussed in the earlier section. At each fixed point there will be an equilibrium 
elongation 5{a) and a corresponding stress S — U*k8{(t) develops in the system (bundle). 
This 5 {a) can be easily expressed in terms of U*{a). This requires the evaluation of c*, the 
internal stress per fiber developed at the fixed point, corresponding to the initial (external) 
stress (7 (= F/N) per fiber applied on the bundle when all the fibers were intact. From the 
first part of eqn. (19), one then gets (for a > ul) 

1-CT* 



= 1 - 



(21) 



(l-(7i) l-GL 

Consequently, 

k5{(t) = (7* = 1 - W{1 - (7l). (22) 

It may be noted that the internal stress a* is universally equal to 1/2 (independent of ctl) 
at the failure point cr = Uc of the bundle. This finally gives the stress-strain relation for the 
RFB model : 



< (7 < ' 
k5{1 — k5) / {1 — a^), (TL<a<ac >■ 

0, (7 > (7c 



(23) 
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Fig. 5: Schematic stress (5')-strain (i5) curve of the bundle (shown by the soUd Une), following eqn. (23), 
with the fiber strength distribution (18) (as shown in Fig. 4). 

This stress-strain relation is schematically shown in Fig. 5, where the initial linear region has 
slope K (the force constant of each fiber). This Hooke's region for stress S continues up to 
the strain value 5 — (Tl/ until which no fibers break {U*{(t) = 1). After this, nonlinearity 
appears due to the failure of a few of the fibers and the consequent decrease of U*{a) (from 
unity). It finally drops to zero discontinuously by an amount a*U*{ac) — 1/[4(1 — (7^)] = a"c 
at the breaking point a = ac oi 5 — (t*/k = 1/2k for the bundle. This indicates that the 
stress drop at the final failure point of the bundle is related to the extent (a^) of the linear 
region of the stress-strain curve of the same bundle. 

Here, the plasticity (nonlinearity) in the response of the bundle comes naturally from 
partial failure of the fibers (and the consequent redistribution of stress among the surviving 
fibers), after the assumed linear region until the lower threshold of failure (18). The 
total failure of the bundle is again discontinuous here and the entire nonlinear response 
characteristics is analytically calculable in this simple model. 
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VIII. Summary and conclusions 

We have reported here an analytic study of the failure dynamics and the consequent plastic 
deformation characteristics of the random fiber bundle model having the property of global 
load sharing. This has been done here for uniform strength distribution p{a) of fibers in the 
bundle (up to a cutoff). As mentioned before, this has been possible due to the inherent 
mean-field nature of the model. The recurrence relation (1) captures essentially all the 
intriguing features of the dynamics. We found that both the breakdown susceptibility x and 
the relaxation time r diverge as the applied stress a approaches its global failure point cXc 
(=1/4 for the uniform strength distribution as shown in Fig. 2) from below, with the same 
exponent value 7 = a = 1/2. The critical dynamics of the model follows a strict power 
law decay at cr = (7^. Ut — 1/2 oc t~^. Though we have identified O = U*{(t) — 1/2 as the 
order parameter (with exponent (3 = 1/2) for the continuous transition in the model, unlike 
conventional phase transitions it does not have a real-valued existence for a > ac- From 
finite-size scaling study, we see that there is a correlation length which diverges with an 
exponent = 1, as (Xc is approached from below. The avalanche size distribution P{m) for 
this mean-field dynamics of the RFB model is given by P{m) ~ , t] = 3. This has been 
confirmed here numerically. As mentioned before, this result is valid for the avalanche sizes 
defined through (16), where the external load on the bundle increases uniformly until the 
total failure at ac- The present as well as the earlier numerical results H, |^, all confirm 
that the analytic results for the exponents a, 7 and r] (for r, x and P(m) respectively) 
are not necessarily restricted to the uniform distribution of fiber strength (assumed here) 
and are more generally valid. The model also shows realistic plastic deformation behavior 
with a shifted (by ctl, away from the origin) uniform distribution of fiber strengths. The 
stress-strain curve for the model clearly shows three different regions: elastic or linear part 
(Hooke's region) when none of the fibers break {U*{a) = 1), plastic or nonlinear part due 
to the successive failure of the fibers {U*{(r) < 1) and then finally the stress drops suddenly 
(due to the discontinuous drop in the fraction of surviving fibers from U*{ac) to zero) at 
the bundle failure point at a = ac (= 1/[4(1 — ai)] for the failure strength distribution 
(18)). Simplicity of the model and consequently of the recurrence relation for the breaking 
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dynamics allows it to have exact analytic results for all its static and dynamic behaviors of 
breaking and the resulting plasticity. 
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